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We study Fermi liquid instabilities in spin-orbit-coupled metals with inversion symmetry. By in¬ 
troducing a canonical basis for the doubly degenerate Bloch bands in momentum space, we derive the 
general form of Landau interaction functions. A variety of time-reversal-invariant, parity-breaking 
phases is found, whose Fermi surface is spontaneously deformed and spin-split. In terms of sym¬ 
metry, these phases possess gyrotropic, ferroelectric and multipolar orders. The ferroelectric and 
multipolar phases are accompanied by structural distortions, from which the electronic orders can be 
identified. The gyrotropic phase exhibits a unique nonlinear optical property. We identify correlated 
electron materials that exhibit these parity-breaking phases, including LiOsOa and Cd2Re207. 
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Novel physics from strong spin-orbit coupling in quan¬ 
tum materials is currently attracting widespread interest 
across many disciplines in condensed matter physics. In 
particular, there is now an intensive investigation of the 
interplay between spin-orbit coupling and electron corre¬ 
lation in d-orbital and /-orbital systems[l|-@]. The ma¬ 
jority of studies have been focused on correlated band or 
Mott insulators, whereas spin-orbit coupling in correlated 
metals has received less attention. It is well-known that 
spin-orbit coupling in metals without inversion symme¬ 
try generates spin-split energy bands and spin-polarized 
Fermi surfaces [4J. This has interesting consequences in 
the presence of electron-electron interactions [5j [l 1]. In 
contrast, in metals with inversion and time-reversal sym¬ 
metry, Bloch bands are doubly degenerate everywhere in 
momentum space. The effect of spin-orbit coupling is 
more subtle: it leads to spin-orbit-entangled Bloch wave- 
functions, which have different spin-polarizations on dif¬ 
ferent atomic orbitalsjri, 3). For this reason, the impor¬ 
tance of spin-orbit coupling in inversion-symmetric ma¬ 
terials can be easily overlooked. 

In this Letter, we explore the consequences of having 
both strong spin-orbit coupling and electron interaction 
in metals with inversion symmetry. By generalizing Lan¬ 
dau’s Fermi liquid theory to spin-orbit-coupled metals, 
we theoretically predict a variety of new ordered phases 
resultin g fr om Pomeranchuk type instabilities in the spin 
channel [14[, which spontaneously break inversion sym¬ 
metry. These phases can be regarded as new examples of 
electronic liquid crystals [l5j , which preserve the trans¬ 
lational invariance and break the point group symme¬ 
try of the lattice. Importantly, because of the spin-orbit 
coupling, these phases exhibit spin-split Fermi surfaces 
with characteristic spin textures, and the onset of elec¬ 
tronic parity-breaking orders is generally accompanied by 
structural changes. We focus on three different parity¬ 
breaking phases having the symmetry of ferroelectric, a 
multipolar and an isotropic gyrotropic liquid respectively, 
and identify their realizations in correlated electron ma¬ 
terials. 


Landau’s Fermi liquid theory of metals starts from 
Bloch states on the Fermi surface. In the presence of 
spin-orbit coupling, Bloch states are not spin eigenstates, 
but remain doubly degenerate at every k in systems with 
both time-reversal ( T ) and inversion (P) svmmetrv(l^|. 
To develop Fermi liquid theory of such spin-orbit-coupled 
systems, we must first choose a basis {|'0k,i), IV^)} for 
the degenerate bands over the entire Fermi surface. As 
observed by Blount long ago[16|, due the absence of spin 
conservation, the choice of basis is not unique: an arbi¬ 
trary U (2) rotation on the doublet at every k produces a 
new basis that appears to be as good as the old one. This 
leads to significant complications, as the form of Landau 
energy functional is basis dependent. 

In this work, we introduce a canonical basis that we call 
“manifestly covariant Bloch basis” (MCBB). This basis is 
defined universally and uniquely by demanding the Bloch 
wavefunctions at r = 0—a two-component spinor—to be 
fully spin-polarized along a global spin-quantization axis: 


0k,i(r = O) = Uk | /), 

0k,2(r = O) = Uk| /). (1) 


where u k and Uk are real and positive ; f, 4- labels elec¬ 
tron’s spin. Importantly, the origin of real space coor¬ 
dinate r = 0 is chosen to be the center of point group 
symmetries of the crystal; and the condition HD is im¬ 
posed on Bloch states on the entire Fermi surface. 

The MCBB can be explicitly constructed by start¬ 
ing from an arbitrary basis {|</k,i), | </>k, 2 ) }- Because 
of time-reversal (T) and inversion (P) symmetry, the 
two members form a Kramers doublet under the com¬ 
bined operation PT[lJ], so that with a proper choice of 
phase, we have |0k,2) = PT |</k,i)- Therefore, the cor¬ 
responding spinors defined by Bloch wavefunctions at 
the inversion center form a Kramers doublet under T: 
0k ,2 (r = 0) = T0 ; i(r = 0), and thus are orthogonal. 
This orthogonality condition guarantees one can perform 
a 17(2) transformation on {|0k,i), |0k,2}} to obtain a new 
basis satisfying CD, or equivalently MCBB. 
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The advantage of MCBB lies in its remarkably simple 
transformation property under point group symmetries, 
which act on both the electron’s spatial coordinate and 
spin. For a generic choice of basis, a symmetry action 
G will map Bloch states at k into those at Gk (or the 
star of k) up to a complicated, k-dependent (7(2) basis 
transformation [ 1 fij . In contrast, the defining property 0 
guarantees that G maps the MCBB IV'k.a) at k directly 
to its partner at Gk, 

G : \tpk, a ) —t U a p(G)\il>Gk,p) (2) 

where (7(G) is the 5(7(2) matrix representation of G. 
Furthermore, the MCBB at ±k are related by time re¬ 
versal symmetry in the same way as spin eigenstates: 

T\ifoi, a ) = e a ,p\^-k,p}- (3) 

Eq. 0 and 0 show that the two members of the MCBB 
a = 1,2 transform identically as spin up and down un¬ 
der symmetry operations. Therefore, for the simplicity 
of presentation, we will refer to the a index of MCBB 
as spin, with the understanding that |^k,a) are not spin 
eigenstates. MCBB provides the starting point for our 
Fermi liquid theory, and is expected to have wide appli¬ 
cations in spin-orbit-coupled systems in general. 

Fermi liquid theory relates the change of energy SE to 
the change in the distribution function of Bloch quasi¬ 
particles up to second order. The distribution function is 
a 2 x 2 Hermitian matrix in spin space, which we write as 
n a p{ k) in MCBB, i.e., n a p{ k) = (cj^Ck,^)- SE is then 
a quadratic functional of n a| a(k), where k is near the 
Fermi surface. For spin-orbit-coupled systems, we find it 
convenient to decompose n a p{ k) in terms of the density 
and spin distribution function: 

^a/3<k) = n(k)6 a p + s(k) ■ <j a p- (4) 

Based on symmetry considerations, we now relate SE to 
the change in the density and spin distribution function. 
First, note the transformation property of n(k) and s(k) 
under time reversal and inversion, 

T : n(k) — > n{— k), s(k) —» — s(—k) 

P : n(k) —¥ n(—k), s(k) —» s(—k). (5) 

It follows that when both symmetries are present, SE 
consists of density-density interaction and spin-spin in¬ 
teraction, taking the form of 

SE = ekdn(k) + F n ( k, k')dn(k)dn(k / ) 

k k,k' 

+ ]Tf-( k,k') Si (k) S ,(k'). (6) 

k.k' 

The invariance of SE under crystal symmetry transfor¬ 
mations further constrains the momentum dependence of 


the interaction functions F n (k, k') and F’ s (k, k'). Un¬ 
like spin-rotationally-invariant systems where the spin 
interaction is isotropic in spin space (7)* oc d l: j ), both 
7 s (k, k') and 7*-(k, k') in spin-orbit-coupled systems are 
constrained by crystal symmetries acting on electron’s 
coordinate and spin in combination. It follows from the 
symmetry property of MCBB © that under a crystal 
symmetry operation G, n(k) and s(k) transform as a 
scalar and a vector field respectively 

G : n( k) ->■ n(Gk) 

s i (k)-^G ij -« i (Gk). (7) 

where G^ is the 50(3) matrix representation of G. 
Hence, F n { k, k') and F^( k, k') transform as a scalar field 
and a rank-two tensor field respectively: 

G: F n (k,k') 7 n (Gk,Gk') 

^(k,k') Gu'Gjj'F^j, (Gk,Gk'). (8) 

Eq.© and © give the general form of the energy func¬ 
tional of Fermi liquids in spin-orbit-coupled systems. 
The novelty here lies in the spin interaction, which is 
anisotropic in spin space and spin-momentum locked. We 
will now explore consequences of spin interactions in spin- 
orbit-coupled Fermi liquids. 

To proceed, we write the spin interaction in a separable 
form given by products of basis functions of k and of k': 

5E spin = k,k>i(k) aj (k') 

k.k' 

= EE F r) ^(k,s(k))<() r) (k , ,s(k')). (9) 

r? k.k' 

Naturally, different basis functions ^(k) fall into differ¬ 
ent representations of crystal symmetry group. As a first 
step, it is instructive to start from isotropic spin-orbit- 
coupled liquids with the largest symmetry group 50(3), 
invariant under any arbitrary rotation of space and spin 
taken in combination. Then, the basis functions are la¬ 
beled by three quantum numbers r) = L is 

the orbital angular momentum, J is the total angular 
momentum J = L + S with 5=1, and J z is the z- 
component J. The 27 + 1 basis functions with the same 
(7, J) and J z = — J, — J + 1,... J form a multiplet. As¬ 
sociated with each (7, J)-multiplet is an interaction pa¬ 
rameter F[ j , which parameterizes the excitation energy 
of a particular type of Fermi surface deformation. Differ¬ 
ent (7, J)-multiplets correspond to “orthogonal” modes 
of Fermi surface deformations. 

We now explicitly decompose the spin interaction into 
a few lowest (7, J)-multiplets. Up to 7 = 1, there are 
four multiplets: (7 = 0; J = 1) and (7 = 1; J = 0,1, 2), 
and hence <5T? sp i n takes the form: 
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SE s pin = ]T Fq s(k) • s(k') + F ? (k • s(k)) (k' • s(k')) + (k x s(k)) • (k' x s(k')) + F 3 S Q tf (k)Qy(k'), (10) 

k,k' 


where Q ,; 7 = Qji is a second-rank tensor constructed from 
k and s(k): 

Qij( k ) = \ (^Sj(k) + k jSi ( k)) - ik • s(k )Sij. ( 11 ) 

Eq. (flOl) is a main result of this work, which shows the 
presence of three p-wave (L = 1) spin interaction chan¬ 
nels parameterized by F(, F£ and F 3 . Similar decomposi¬ 
tions of spin interaction into higher angular-momentum 
channels can be carried out by constructing high-rank 
tensors from powers of k and s(k). 

At this point, it is worth discussing the effect of peri¬ 
odic crystal potential, which reduces the full rotational 
symmetry to its subgroup, the point group of a crys¬ 
tal. In this case, Fermi liquid interactions can still be 
decomposed into different spin-orbit-coupled channels as 
in m- However, these channels are in one-to-one cor¬ 
respondence with the irreducible representations of the 
point group, instead of the (L, J)-multiplets for S , 0(3) 
group. Despite this difference, for many crystal struc¬ 
tures such as cubic, tetragonal, trigonal and hexago¬ 
nal, the four channels in m remain to be in differ¬ 
ent point group representations, and hence orthogonal 
to each other. 

When one or more interaction parameters in the spin 
channel become negative and of sufficiently large mag¬ 
nitude, Fermi surface instability occurs. A well-known 
example is the ferromagnetic instability associated with 
Fq in the s-wave spin channel. This work is concerned 
with Fermi liquid instabilities in the p -wave (more gen¬ 
erally, odd L) interaction channels in spin-orbit-coupled 
metals. The resulting phases are time-reversal-invariant 
and parity-breaking, and as we will show, exhibit novel 
properties arising from spin-orbit coupling. 

First, consider the instability associated with F* in 
the (L = 1; J = 0) channel. According to (flTil) . this 
instability generates an Ising order parameter: 

7? = ^k-s( k ). ( 12 ) 

k 

r] is a pseudo-scalar, because it is invariant under time- 
reversal and all rotations, but breaks inversion and all 
reflections. Therefore the ordered phase with rj ^ 0 is 
an isotropic gyrotropic liquid. This gyrotropic order pa¬ 
rameter splits the original spin-degenerate Fermi surfaces 
into two with unequal volumes, with opposite spin polar¬ 
izations. Unlike the case of ferromagnetism, here the spin 
quantization axis defined in terms of MCBB is not uni¬ 
form but parallel to the momentum: s(k) oc 77 k, which 
leads to a hedgehog spin texture over the Fermi surface. 


Next, consider the instability associated with F 2 S in 
the (L = 1; J = 1) channel. According to Ea. (flTil) . this 
instability generates a vector order parameter 

P = ]Tkxs(k). (13) 

k 


We observe that P has the same symmetry as the ferro¬ 
electric polarization: it is odd under inversion and invari¬ 
ant under time-reversal, and transforms as a vector under 
rotation. Therefore we identify the ordered phase with 
P ^ 0 as a “ferroelectric” metal that spontaneously de¬ 
velops a polar axis, despite that its charge polarization is 
screened by free carriers[17]. In this phase, Fermi surfaces 
are spin-split and deformed by a spontaneously generated 
spin-orbit field h(k) acting on the original Fermi surface. 
h(k) is k-dependent and proportional to the spin po¬ 
larization field s(k) generated by the ferroelectric vector 
order parameter: h(k) oc s(k) oc P x k. This spin-orbit 
field has the same form as the Rashba spin-splitting due 
to an external electric field. In our case, both the Rashba 
spin splitting of the Fermi surface and the accompany¬ 
ing ferroelectric order P are caused by strong electron 
interactions. 

Lastly, consider the instability associated with F£ in 
the (L = 1, J = 2) channel. According to ([Tfl]) . the corre¬ 
sponding order parameter is a traceless symmetric matrix 
given by 


Qjj — Qij(k)) (14) 

k 

where Qij{ k) is defined in (fill) . This order param¬ 
eter has d -wave symmetry and odd-parity, and hence 
is a second-rank pseudo-tensor. The multipolar phase 
with Qij ^ 0 can be regarded as an electronic analog 
of the chiral nematic liquid crystals fl8| . Its Fermi sur¬ 
faces are spin-split and deformed by the spin-orbit field 
hi( k) oc Si(k) oc Qijkj. If the matrix may have 
two degenerate eigenvalues, the ordered phase is uniax¬ 
ial; otherwise, it is biaxial. 

To summarize, we find new parity-breaking phases 
with gyrotropic, ferroelectric and multipolar orders in 
spin-orbit-coupled Fermi liquids, driven by strong p -wave 
spin interaction. Because of spin-orbit coupling, these 
symmetry-breaking phases exhibit spin-spilt Fermi sur¬ 
faces. The magnitude and direction of the spin-splitting 
vary strongly over the Fermi surface. The characteris¬ 
tic Fermi surface splitting and spin texture in momen¬ 
tum space are predicted to be the hallmark of parity¬ 
breaking phases in spin-orbit-coupled metals. These fea- 
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tures of parity-breaking phases can be detected by angle- 
and spin-resolved photoemission spectroscopy. 

Since the above parity-breaking order parameters are 
time reversal invariant and break rotational symmetry of 
the crystal, they couple linearly to lattice distortions that 
lowers the point group to the same subgroup without en¬ 
larging the unit cell. As a result, the transition driven 
by Fermi liquid instability is generally accompanied by 
a structural transition, from which the electronic order 
can be inferred. A possible exception is the gyrotropic 
order m, which preserves the full rotational symme¬ 
try of the crystal. Such a high degree of symmetry may 
not be compatible with any lattice distortion caused by 
atomic displacements. In this case, it will be difficult to 
detect the electronic gyrotropic order with conventional 
methods 25-281]. 

It is important to emphasize that the presence of spin- 
orbit coupling is indispensable to the parity-breaking 
phases found in this work. To make this point clear, let us 
consider spin-rotationally-invariant Fermi liquids, whose 
spin interaction in the p -wave channel takes the form 
^ kk , (k • k') s(k) • s(k'). The corresponding Fermi liq¬ 
uid instabilities have been studied in detail 1§- 211. The 
ordered phases were found to simultaneously break two 
symmetries, the rotational symmetry of space and of 
spin. As a result, the spin textures are free to rotate 
as a whole, instead of being rigidly locked to momentum 
as in our case. Moreover, these electronic orders cannot 
couple directly to lattice distortions (which preserve spin 
rotational symmetry), unlike the ferroelectric and mul¬ 
tipolar phases of spin-orbit-coupled Fermi liquids. We 
also note that besides spin-orbit coupling, dipolar in¬ 
teractions in ultracold Fermi gases, which also lock spin 
and momentum, can generate ordered phases with simi¬ 
lar features 22|, [23[. In addition, parity-breaking phases 
can occur in spin-orbit-coupled insulators [23]. 

Finally, based on recent experiments, we identify sev¬ 
eral correlated electron materials that show evidence 
of the above parity-breaking orders. First, a recently 
synthesized material LiOsOa was found to undergo a 
second-order ferroelectric structural transition at low 
temperature[29;]. The high-temperature structure is Z? 3 d, 
which is inversion symmetric. A polar axis in c direc¬ 
tion appears in the low temperature phase, reducing the 
crystal symmetry to C^ v . Based on the observation of 
unusually large residual resistivity and Curie-Weiss be¬ 
havior of spin susceptibility, it has been suggested that 
electron correlation plays an important role and possi¬ 
bly drives the structural transition}^. Therefore the 
low-temperature phase of LiOsOa may be an electronic- 
driven ferroelectric metal. 

Second, we suggest pyrochlore oxides A 2 B 2 O 7 as 
promising candidates for the multipolar phase. The py¬ 
rochlore crystal structure has the Oh point group symme¬ 
try. Due to this crystal anisotropy, the five-component 
multipolar order parameter Qij defined in ( 1111 ) splits into 


a two-dimensional E u representation ( Q xx — Q yy , 2 Q zz — 
Qxx ~ Qyy) and a three-dimensional T 2 U representa¬ 
tion (Q xyi Q yz ,Q zx ). Recently, the pyrochlore oxide 
Cd 2 Re 2 07 was found to undergo a second-order struc¬ 
tural transition from cubic to tetragonal at T c = 200K, 
with an order parameter of the E u symmetry [3p| - [33j . Re¬ 
markably, the lattice change across the transition is ex¬ 
tremely small, whereas electrical properties change dras¬ 
tically. In addition, a large mass enhancement above the 
transition temperature was inferred from transport and 
optical measurements [13-0]. Therefore, the structural 
transition in Cd 2 Re 2 C >7 may be induced by an electronic 
transition to the mulitpolar phase. 

In the above examples of ferroelectric and multipolar 
phases, the appearance of electronic order is inferred, by 
symmetry consideration, from the structural distortion it 
couples to. It is desirable to directly probe the change 
in electronic structure, Fermi surface, and spin texture 
across the parity-breaking phase transition via, for exam¬ 
ple, angle- and spin-resolved photoemission spectroscopy. 
Moreover, it will be interesting to determine whether the 
driving force for the transition is structural or electronic. 
On the other hand, the isotropic gyrotropic order can be 
more elusive. In the case of pyrochlore crystals, the gy¬ 
rotropic order belongs to the A lu representation of the 
Oh point group, which is incompatible with any phonon 
mode at the Brillouin zone center 371]. Therefore, it can¬ 
not be generated by structural distortions, and if found, 
should have an electronic origin. 

In addition to photoemission, nonlinear optics is a 
powerful tool for detecting parity-breaking orders de¬ 
scribed in this work. For example, the multipolar phase 
in tetragonal Cd 2 Re 2 C >7 has been successfully detected 
by second-harmonic generation (SHG)0|. Regarding 
the gyrotropic order (1121) . we find it has the same sym¬ 
metry as the rank-3 isotropic tensor e^fe. It then follows 
from symmetry that this gyrotro pic order should lead to 
sum-frequency generation (SFG)[39(, in which two inci¬ 
dent fields Ej _2 at different frequencies generate an 
electric dipole P at the frequency u> = u>i + u> 2 ' 


Piiyj) oc sgn.(T))eijkE lt j(uii)E2,k(u2)- (15) 

This nonlinear optical effect gives a direct way of detect¬ 
ing the much-hidden gyrotropic order. 

Our work on spin-orbit coupled metals leaves a number 
of open questions for future studies. It is worthwhile to 
relate the phenomenological parameters in Fermi liquid 
theory to microscopic interactions. It will be extremely 
interesting to study superconducting instabilities of spin- 
orbit-coupled Fermi liquids, especially those proximate to 
the parity-breaking phases. 

Note: A recent work shows that optical circular dichro- 
ism can be used to probe the parity-breaking phases de¬ 
scribed in this work[40j. 
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